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Abstract 

Hyperspherical partial wave theory has been applied here in a new way in the cal- 
culation of the triple differential cross sections for the ionization of hydrogen atoms 
by electron impact at low energies for various equal-energy-sharing kinematic condi- 
tions. The agreement of the cross section results with the recent absolute measure- 
ments of Roder et al [51] and with the latest theoretical results of the ECS and CCC 
calculations [29] for different kinematic conditions at 17.6 eV is very encouraging. 
The other calculated results, for relatively higher energies, are also generally satis- 
factory, particularly for large 60b geometries. In view of the present results, together 
with the fact that it is capable of describing unequal-energy-sharing kinematics [35], 
it may be said that the hyperspherical partial wave theory is quite appropriate for 
the description of ionization events of electron- hydrogen type systems. It is also clear 
that the present approach in the implementation of the hyperspherical partial wave 
theory is very appropriate. 
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I. INTRODUCTION 



In the study of electron hydrogen atom ionization cohision, the simplest three-body 
ionization problem in atomic physics, there are many attempts for a complete solution 
but all of these face tremendous difficulties and have only limited success. Except for one 
or two attempts all use time-independent framework. For accurate information regarding 
scattering events, one may solve accurately the Schrodinger equation for the scattering 
states '^l'^^ or [see Newton [1] for their definition] given by 

H^iy = My (1) 

taking account of the appropriate boundary conditions. 

Ionization amplitudes may then be obtained either from the flux condition at infinity 
or from appropriate projections. In the literature both '^[^^ and ^'^ ■* have been widely 
used. There are a large number of attempts which strives to solve for ^'^-^^ Among these 
the most successful attempts are the various close-coupling calculations [2-4]. In these 
calculations '^[^^ are expanded in terms of basis functions and ionization information 
are extracted from a solution of the unknown expansion functions. Another possibility 
is to expand in terms of a complete set of functions in the angular variables. In 

these regards the attempts of Kato and Watanabe [5, 6] are remarkable. They used 
hyper spherical co-ordinates and expanded ^^"""^ in terms of hyper-radius dependent an- 
gular functions. Matching with a wave function, which satisfies an approximately correct 
boundary condition, they obtained with remarkable success, the total ionization cross 
sections down to the threshold. However, differential cross section results of this the- 
ory are not known. Very recently Rescigno and associates made [7, 8] a breakthrough 
calculation and reproduced for equal energy-sharing and constant angular separation 
Qab of the outgoing electrons, the cross section results, at low energies, with surprising 
success. In these calculations they expanded 'I'^^'* in terms of spherical harmonics in 
four angular variables. Then they converted the resultant differential equations for the 
radial functions, in two radial variables, into a set of difference equations over a large 
network in the radial variables-plane. They used a novel technique. Using a complex 
scaling procedure they converted the scattering problem as if into a bound state prob- 
lem. Then they solved a huge set (several million) of linear equations using very special 
techniques. Ultimately they obtained ionization amplitudes using the flux condition. 
Later [9] they confirmed their results using projection technique. Although the ECS ap- 
proach reproduced the equal-energy-sharing, constant-Oa;, results perfectly well, results 
of this approach for unequal-energy-sharing kinematics are not known. There are also 
large number of attempts of using ^ in extracting ionization information. In such 
cases projection approach has been generally used. There the ionization amplitudes are 
calculated from 

T!^ = {^%'^\V^\^^). (2) 
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Brauner, Briggs and Klar [10] and later Berakdar [11] and Berakdar et al [12, 13] made 
use of this approach. They used ^'^"^ which are asymptotically correct (or nearly so) 
but are unlikely to be correct at finite distances. As a consequence results of these 
calculations are only moderately accurate. Moreover there are no systematic tractable 
way of improving the results. 

An alternative approach for determining the electron atom collision cross sections 
is to solve a coupled set of integral equations for the off-shell T-matrix elements. Das 
and associates [14-16] have used this approach in the study of various electron hydrogen 
atom and electron helium atom scattering problems by solving the resultant equations 
in a rather crude manner. However, they always obtained moderately good results. 
There are also attempts [17-24] to improve the calculations. Along these lines the most 
successful calculations are the convergent close coupling (CCC) calculations of Bray et al 
[25-30]. In many contexts they applied the CCC method with surprising success. Their 
latest results [29, 30] claim accuracies comparable with the ECS results at low energies. 

Another promising approach for the electron hydrogen atom ionization problem is 
the hyperspherical partial wave approach [33-34]. Details of this approach are given in 
[34]. In section II we also present important features of this approach. Earlier with an 
additional approximation of neglecting the coupling effects, some results were obtained 
[35, 36] which are qualitatively not very bad. Recently this approach has been used 
[37, 38] retaining fully the coupling effects. In solving the relevant coupled set of radial 
wave equations over an initial interval [0, A], R-matrix [39] approach had been used. 
Although the results were always found to be of the correct magnitude, pseudo-resonance 
type behaviour gave much troubles in extracting correct cross section results. To avoid 
this problem we use a new approach. This appears to be very successful and leads to 
very interesting results both for equal-energy-sharing constant-0ab kinematics, equal- 
energy-sharing asymmetric kinematics, and also for unequal-energy-sharing kinematics 
[38]. Thus it appears that hyperspherical partial wave theory is quite appropriate for the 
study of ionization problems of electron-hydrogen type systems. 

Most recently two very broad-based theories have been proposed. One of these is the 
time-dependent close coupling theory [40] and the other is the hyperspherical R-matrix 
theory [41]. Positions of these theories are not yet very clear. 

II. HYPERSPHERICAL PARTIAL WAVE THEORY 

In the hyperspherical partial wave theory one uses the time-independent framework. 
In the time-independent framework the T-matrix element is given by expression (2) or 
alternatively by 

Tf, = i^flvM-^^). (3) 

In these expressions <I> j and <I> j are the unperturbed initial and final channel wave func- 
tions, satisfying certain exact boundary condition at infinity and that Vi and Vf are 



3 



the corresponding perturbation potentials. For the case of ionization of hydrogen atoms 
expression (2) is more appropriate for use, since in this case asymptotically correct $j 
is easily available. Many use expression (3), including ECS [9] by projection method, 
but inappropriately, since the corresponding ^/'s they use do not satisfy the correct 
boundary condition. In the hyperspherical partial wave theory ^'^ is expanded in 
terms of hyperspherical harmonics, which arc functions of five angular variables. The 
corresponding radial waves are functions of one radial variable, the hyperradius R only. 
This proves to be advantageous in numerical computations, since then the five angular 
variables range over a bounded compact domain, while only one variable R ranges over a 
semi-infinite domain [0, oo). It may be noted here that so far nobody could take account 
of the exact boundary condition in the asymptotic domain for the accurate solution of 
^j"-*. Here we aspire to take account the exact boundary condition at infinity, in the 
limit. This is the most novel feature in the hyperspherical partial wave theory. Here we 
may note that two plane waves exp{ipa • ri)/(27r)^/^ and exp{ipi, ■ r2)/(27r)^/^ may be 
decomposed in partial waves as usual and then these may be combined (using a formula 
in Erdelyi [42]) to obtain an expansion in terms of hyperspherical harmonics (j)\{u)), in 
five angular variables w = (a, ^i, ^i, 02)- A symmetrized two-particle plane wave has 
the expansion [Das, 34] 

[exp{ipa ■ rl + ipb ■ f^) + {-iyexp{ipb ■ fl + ipa ■ r2)]/(27r)^ 



where ux = X + ^ and X = li + I2 + 2n (A also denotes the multiplet (/i, l2,n) depending 
on the context). Here R = \/rf + r|, a = atan{r2/ri), fl = (ri,6i,(pi), = {r2,02,4'2)- 



Similarly P = ^pI+pI, oq = atan{pb/pa), Pa = {Pa,da,(pa), Pb = {Pb,Gb,4>b), and 
p = PR, and = ("0, Oa, 4>a, Ob-, 4>b), 

= i{l + (-l)-^+^+"}P,7(a)3^,f^(ri,r2)}, for h = h = I, (5) 

and a corresponding expression for (pxi^o) (similar expressions may be easily derived for 
product of more than two plane waves). 

Now the symmetrized wave '^^fj may be expanded in terms of symmetrized hyper- 
spherical harmonics (?!)^'s as 

*(;)(i?,u;) = 2^1^ ^ cPiiu), (6) 

(s) 

where satisfy an infinite coupled set of equations 



Ft\R)+j:^-^Fif\R)=0, (7) 
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for each symmetry s (s = for singlet and s = 1 for triplet)and for each total angular 
momentum L (and its projection M, and so also for a definite parity tt). In the above 
expression 

ciy = -{•H\C\<l>'y)/P, and 



+ 

cosa sma \ricosa — r2Sina\ 
The equations (7) have to be solved over the semi-infinite domain [0, oo). Solution in 
the asymptotic domain is simple. The equations have nice asymptotic solutions. One 
may note that because of conservation rules the equations (7) are decoupled into sets for 
fixed ijl = {L,S,tt) and different N = (Zi,Z2,n). So we set F^^^ = and, for the set 
with fixed /x (and dropping fi from f^) we can write equations (7) as 



dB? 



JN + 7^ JN' 



N' 



R 



0, 



for a solution over the finite domain and 

l^N {l^N + 1) 



f. 



N' 



0, 



(8) 



(9) 



for solution over an asymptotic domain, say [Ro,oo). Next we consider the solution 
problem first over an asymptotic domain [Rq, oo) and then over the finite domain [0, Rq]. 



A. Solution in an asymptotic domain 

The equations (9) have two sets of solutions [34] in an asymptotic domain [i?o,oo) of 
the form given by 

oo (^) . ^ oo ,(£) n 

/.aw = E^^^+E^^. (10) 

oo (£) . ^ oo m a 

&iP) = E "-"^ + E (11) 

where 6^ = p + a^ln 2p and is the k-th eigen value of the charge matrix A = {a^N') 
and that the coefficients a[?^, c^']^ and df^^ are determined from recurrence rela- 

tions. Thus the coefficients afj^ and bfj^ are determined from the relations - 

2[{Akf + /2/]xJ^') = [AkAk - l{l - l)Ak - l{2l - l)afe/]xf-') 

•2(1 _ nnvC^-i) 



and 



-[(2/ - l)akAk + /Afe - l\l - 1)/]Y^'-^^ (12) 



2[{Akf + i'l]Yl> = [Akhk-l{l-l)Ak-l{2l-l)akI]Yl~'-^ 



+ [(2l - l)akAk + lAk - l\l - m'^k (13) 
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where the N-th components of vectors X^'^ and Y^'^ are given by 
and where 

Ak = A- aki, {A.k)NN' = [ttfe^ + + '^)]Snn'- 

The initial vectors X^°^ and Y^^ are given by 

-^k - -^k, » jfc - U, 

Xjfc being the k-th eigen vector of the charge matrix A and I is the identity matrix. 
Solution for c^'^'s and rfj?]^'s are similarly obtained from the above recurrence relations 
after setting X^°^ = and Y^^ = Xj.. In this way we get solution vectors f^^]^ and f^^^ 
of equations (10) and (11). 

B. Solution over a finite domain 



Here we consider the solution of equations (8) over a finite domain [0, Rq] . Away from 
the origin, solution of the equations is easy. A Taylors series expansion method suffices 
for getting arbitrarily accurate solutions. The main difficulty is in starting the solution 
from the origin. Near origin the equations have analytic solutions ( Fock [43] )but these 
are too complicated to be useful in numerical computations. In our earlier calculations 
[37, 38] we used R-matrix approach for getting solutions over an initial interval [0, A] 
(with A suitably chosen). But this approach faces difficulties as pseudo-resonance type 
behavior appears giving much troubles in determining appropriate solutions. To avoid 
such difficulties we consider here a new approach. For the interval [0, A] we consider a 
boundary value problem . The solution vector fQ''^(R) is assumed to have a value at 
the origin and the k-th column of the N„ix x N„ix identity matrix at i? = A. We divide 
the interval [0, A] into m subintervals and use a five-point difference formula 

fN{Rk) = Y^[-fN{Rk-2) + WN{Rk-l) - 30/jv(i?fe) + 16fN{Rk+l) - fN{Rk+2)] 

+{P^n\0} (14) 



for A; = 2, 3, ■ ■ ■ , m — 2, and a formula 

fN{R + h') = 2f';,{R + 2h')-f';,{R + 3h') + ^[fN{R)-4.fN{R + h') + 6fN{R + 2h') 

-4MR + 3h') + MR + 4/i')] + {-^/('"■^(O}- (15) 

with R = Rq, h' = h ioT the equation sd, R = Ri and R = Rm , h' = —h ioi the equation 
at i? = i^TO-i.For continuing solutions from A onward we need first order derivatives at 
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A. For this we use the simple formula 



fNiRm) = [-fN{Rn.-4)+24fN{R„,-2)-128fN{Rm-l) + W5fN{Rm)]/{Mh) 

+ ffN{Rm) + {-^f^^\0}- (16) 

The resultant linear set of equations ultimately take the form 

SYC^) = b, (17) 

where is a sparse matrix. We solve these equations by the biconjugate gradient method 
[44] using routines given in [45] . With a suitable choice of a preconditioner the iterations 
smoothly converge (with a few hundred iterations or even less) to five or six decimal 
places for a suitable choice of error limit (say, 1 in 10^ parts). In this way N^x solution 
vectors fg"*^ are determined over [0, A]. The solutions are next continued over [A, i?o] 
by Taylors expansion method with stabilization [47] after suitable steps, giving solution 
vectors over [0, i?o]- 

C. Matching of the solutions: Determination of ^fy, ^ 

For finding the physical solution vectors fph and the scattering state we first 
define solution matrices fo, fgn and fcs, by putting side by side the corresponding solution 
vectors fp^^ fin\ fcs\ for A; = 1, 2, • • • , N^n^- Then the physical solution vector fph may 
be defined over [0, Rq] by 



k=l 



wr)=j: 9kfS'\R) (18) 

and over [Ro,oo) by 



fph(i?) = ^''"''C,4n\p) + ^'^'-^d.fis^lp). (19) 
k=l k=l 

2Nmx of the SNjjix unknown coefficients are now determined by matching values (and 
first order derivatives) of the two sets of solutions at a point Rq where all the solutions 
are valid. The remaining N^x unknown coefficients are then determined from the de- 
mand that "^^fj actually satisfies the appropriate boundary condition. To facilitate the 
computations we first define the K-matrix through the relation 

h-B = fsn + fcs ■ K (20) 

where B is some unknown constant matrix. (The K-matrix thus defined is a little different 
from the one usually defined. But in any case it should be symmetric.) The K-matrix 
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is then determined by matching values and first order derivatives of the two sides of 
equation (20). Then in the asymptotic domain one has 



where 



fph = {fsn + fcs- K) -C 
~ fsn ' C ~l~ fcs ' d, 

d = Kc. 



(21) 



(22) 



Finahy fph is completely determined once the vector c is determined. Now c is deter- 
mined from the consideration that ^'^"^ is asymptotically a distorted plane wave (repre- 
senting the two outgoing electrons) plus incoming waves only. So we equate coefficients 
of the outgoing wave exp{ip) of both "^^fj and the symmetrized plane wave (4) (except 
for the distorting terms exp{iakln2p). This gives 



c = 



I + iK 



-1 



(23) 



where 



where X is the matrix consisting of eigen vectors of the charge matrix A (and is non- 
singular) and ^^^^*{ujq) is given by 



/ 0^*1 (^o) \ 



*^*(a;o) 



V 



N, 



(24) 



In this way the physical radial vectors fph(-R) are determined for each fi = (L,S,7r) and 
ultimately the full (but approximate) scattering state "^^fj is obtained. 

Substituting this expression in equation (2) one obtains the scattering amplitude in 
the form 

r(a;o) = ^T;, = i-^^"'^C^(7V)</,|,(a;o) (25) 



2tt 



N 



The triple differential cross section is then given by 

d^(^ _ (27r)Vpb j-l (o),2 ^ 3 



dEad^ad^b 



Pi 



fi 



+ 



(26) 



By increasing the number of channels for each /j, = {L,S,Tr) one may expect to 
obtain converged cross section results. 



III. PRESENT CALCULATION 

In our present calculation there are two important parameters A and Rq which are needed 
to be judiciously chosen. The parameter A of the initial interval [0, A], for a solution of 
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the radial equations (8), has been chosen to be 5 a.u. for all the energies considered here. 
The results do not depend on A for small variations (of a few a.u.) about this value. On 

the other hand the choice of the parameter Rq, the asymptotic range parameter, is very 
crucial. Without its appropriate choice the asymptotic series solutions (10) and (11) are 
unlikely to converge. Here it is found that for convergent asymptotic scries solutions Rq 
is needed to be such that Rq ~ where E is the energy in the final channel. Thus 
for energies of 30 eV, 25 eV, 19.6 eV, 17.6 eV, and 15.6 eV this range parameter Rq 
may be chosen greater than the values 60 a.u., 70 a.u., 90 a.u., 120 a.u. and 150 a.u. 
respectively. We have chosen Rq around these values in our calculations. However for 
the computation of single differential cross section (SDCS) it is necessary for converged 
results to vary Rq, and extrapolate, as in ECS calculation [8] for Rq oo. Our limited 
computational resources restrict us to take single Rq value for each energy.Moreover for 
arbitrary large Rq unwanted errors are likely to make the results erroneous. So some 
optimum choice of Rq has to be made for each energy with a few trials. In our present 
computations this has been done. Next we consider the choice of L values for inclusion 
in the calculations for different energies. For 15.6 cV energy, values of L upto 5 proved 
sufficient. For 17.6 eV calculations values of L upto 7 are found to be necessary. For 
the other energies considered here, values of L upto 9 have been included. For fixed 
{L,S,tt) the number of channels, the number of states with different (h,l2,n) triplets, 
which have been included, were chosen suitably for fairly converged results. In any case 
for fixed (Zi, Z2) pairs n was varied from to 9. In this way convergence with respect to 
n is obtained. The number of (Zi, Z2) pairs, which have been included, are somewhat less 
for L > 2 compared to those in the ECS calculation . These pairs arc chosen more or 
less in the order as in ECS calculations (a little different from those of hyperspherical 
calculations of Kato and Watanabe [6]). However there could be some better choice. For 
lower energies convergence with respect to the increase in channel size is rather smooth. 
It is not so for relatively higher energies of 25 eV or 30 eV. Nevertheless wc have obtained 
nearly converged results in the cases considered. All the results presented here are more 
or less based on 50 channels calculations. Most of the calculations, reported here, were 
done on Pentium -III PCs. Calculations for 15.6 eV energy could not be done on PCs. 
Thus the results for 15.6 eV and 17.6 eV, which are presented here, have been derived 
from calculations on a SUN server. It may be further added here that for 15.6 eV 
the SDCS results show that for equal energy sharing case the calculated SDCS value 
is about twice the expected value, although the calculated total cross section appears 
correct (sec table I). This is unacceptable. In any case the various cross section results 
for this energy have been multiplied by a factor 0.5 before presentation in the figures. 
For other energies, however, we have nearly the correct SDCS values for equal energy 
sharing situations. Calculation on a larger scale with larger values of Rq, and with more 
precise solution of the equations (17), may decide absolutely the normalization question 
of the measured results of Roder et al [48] for 15.6 eV energy. Cross section results for 
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15.6 eV energy have been included here for the sake of completeness. 



IV. RESULTS 

A. Triple Differential Cross Section for Constant 6a(,Geometry 

The triple differential cross section results for equal-energy-sharing constant Qab ge- 
ometries are presented in figures 1(a) for 15.6 eV energy, in figure 2(a) for 17.6 eV energy 
and in figures 3, 4 and 5 for energies of 19.6 eV, 25 eV and 30 eV. In these figures we have 
presented the theoretical results of CCC calculations [26, 28, 29, 30, 51] and of ECS cal- 
culations [8, 50]. Here we have also included the absolute measured values of Roder et al 
[47, 48] for 15.6 eV and the most recent re-measured (with necessary inter-normalization) 
values of Roder et al [51] for 17.6 eV energy. For other energies the measured results [47] 
are only relative and are normalized as in [9] . Our results are generally comparable with 
the ECS results in magnitude. For 17.6 eV our present results appear most interesting. 
These are even somewhat better compared to the ECS and CCC results for @ab = ISO'^ 
and 180*^. For these values of @ab, the 15.6 eV results also appear good , particularly 
in shapes, but confirmation by larger scale calculation is necessary. The 19.6 eV results 
also appear to be very good. For other energies our results appear less satisfactory in 
comparison with the ECS and CCC results. 

B. Triple Differential Cross Sections for Fixed 9a Geometry 

In figures 1(b) and 2(b) we have compared our results for equal-energy-sharing asym- 
metric geometries with absolute measured values of Roder et al [48] for 15.6 cV and 
Roder et al [51] for 17.6 eV, in which one of the outgoing electrons is observed in a fixed 
direction while that of the other one is varied. In these cases we again compare our re- 
sults with the calculated results of ECS and CCC theories. Here our results also appear 
to be quite good, particularly for 17.6 eV in view of the most recent measurements. For 
9a = —30° at 15.6 eV, the peak position of our calculated curves are little shifted to the 
right. Otherwise all the results of the present calculation appear satisfactory. 

C. Triple Differential Cross Sections for Symmetric Geometry 

In figures 1(c) and 2(c) we have presented TDCS results for symmetric appearance of the 
two outgoing electrons relative to the incident electron direction, for 15.6 cV and 17.6 eV 
incident electron energies, for which there are again absolute measured results [48, 51]. 
For 15.6 cV energy our results agree qualitatively with the experimental results. Here a 
70-channel calculation has been found to be necessary. For 17.6 eV our results do not 
appear very good. For 15.6 eV energy both the ECS and CCC theories underestimate 
the cross section results considerably. For 17.6 eV energy, however, both ECS and CCC 
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theories give good overall representation. 

D. Integrated Cross Sections and the Spin Asymmetry Parameter 

The parabolic fitted curves to our computed single differential cross sections data are 
generally close to the ECS extrapolated curves but our raw data which could be calcu- 
lated, as in ECS flux method, away from the two ends of the energy intervals widely 
differ from ECS or CCC (wherever available) curves. However, the computed total in- 
tegrated cross sections, with suitable extrapolation from these are generally good. Here 
in table I we have presented values of integrated cross sections aj = {ag + 3crt)/4 and 
the spin asymmetry parameter A = (as — CFtj/i'^s + 3crt) where Ug and at are the singlet 
and the triplet cross sections, together with values of ECS theory by flux approach [8] 
and those of CCC theory and the experimental values . The integrated cross sections 
agree with the experimentally measured values of Shah et al [52] within about 20% . The 
spin asymmetry parameter A agrees, however, excellently with the measurements [53, 54]. 

Next we note down the shortcomings and difficulties associated with the present ap- 
proach. The first point to note is that it may not be possible in this approach to get 
reliable cross section results for extreme asymmetry, as in ECS flux approach, for one of 
the outgoing electrons sharing very small energy values compared to the other. 'Contam- 
ination with high Rydberg states', as in ECS calculation [8] gives wrong results for finite 
values of Rq in such cases. Extrapolation to i?o = oo may only lead to reliable results 
in those cases. This may require larger computational resources. Another difficulty to 
be noted is the appearance of a few large eigen-values of the charge matrix for large- 
channel calculations. In such cases computational strategies are needed to be reviewed. 
In our calculations this has occurred in a few cases. In such cases we simply cut-short 
in magnitude these one or two large cigen values to the neighboring ones. However a 
better approach may be necessary to tackle such problems. No other difficulties appear 
worth mentioning. For a fully converged results inclusion of more channels (about 100 
or a little more) may be required with appropriate choice of (li, I2) pairs (say, as in ECS 
calculation) and with further stabilization. However these are subjects of further studies 
requiring more computational resources and time. 

V. CONCLUSIONS 

The results of the present calculation fairly display the capability of the hyperspherical 
partial wave theory in representing results for cqual-cncrgy-sharing kincmatical condi- 
tions at low energies. The new approach that has been used in the implementation of 
the hyperspherical partial wave theory appears very appropriate. The computed cross 
section results are observed to be very satisfactory. If one recalls the capability of the 
theory to describe the ionization collisions for unequal-energy-sharing asymmetric kine- 
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matic conditions (as indicated in [38]) also then the capabihty of the hyperspherical 
partial wave theory towards a complete description of the electron - hydrogen - atom 
ionization problem is well demonstrated. Considering the computational facilities used 
(Pentium - III PCs and a SUN Enterpriser 450 server) success of the present calcula- 
tion is appreciable. For fully converged results better computational facilities may be 
required. The theory may easily be applied in the study of ionization of hydrogen-like 
ions with a little change in the definition of the wave function $j and the interaction 
potential V^. The theory may also be extended for application to the double-ionization 
of helium atom or helium-like ions or to other multi electron ionization processes. 
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TABLE 



Table I. Total integrated ionization cross sections (a.u.) and the spin asymmetry 
parameter. The experimental values of cross sections are those of Shah et al [52] (the 
starred numbers arc from extrapolation). ECS results arc from [8] and the CCC results 
are from [4]. In the experimental results of the asymmetry parameter of Crowe et al [53] 
and Fletcher et al [54] presented here, the numbers with superscript + or — denote the 
available result just a little above or below the energy considered (for the exact energy 
values the corresponding references are to be seen). For 15.6 eV energy,ECS results of 
earlier calculation [8] are not available. So for this energy we present results from [9] and 
indicate it so in the table. 



Ei(eV) 


30 


25 


19.6 


17.6 


15.6 


Total Integrated Cross Sections 


Present: 


2.13 


1.82 


1.14 


0.83 


0.49 


ECS: 


1.79 


1.62 


1.09 


0.80 


0.36 [9] 


CCC: 


1.92 


1.57 


1.01 


0.75 


0.38 


Expt.: 


1.81* 


1.55* 


1.00 


0.74 


0.39 


Spin Asymmetry 


Present: 


0.31 


0.41 


0.47 


0.55 


0.48 


ECS: 


0.42 


0.45 


0.51 


0.51 


0.52 [9] 


CCC: 


0.41 


0.43 


0.50 


0.51 


0.53 


Expt.: 












Crowe 


0.28 


0.39- 


0.42- 


0.47- 


0.50- 


Fletcher 


0.31 


0.41- 


0.40+ 


0.50- 


0.48- 
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Figure Captions 



Figure 1(a). TDCS in coplanar equal-energy-sharing constant angular separation Qab 
of the outgoing electrons for incident electron energy Ei = 15.6 eV vs. ejection angle 9a 
of the slow outgoing electron. Continuous curves, present results ; dashcd-curvcs, ECS 
results [9, 50]; dash-dotted curves, CCC results [26, 48]. The experimental results arc the 
absolute measured values of Roder et al [47, 48]. Present results have been multiplied 
by a factor 0.5 (see text). 

Figure 1(b). TDCS in coplanar equal-energy-sharing geometry for incident electron 
energy = 15.6 eV for fixed 9a and variable 9b of the out going electrons. Continuous 
curves, present results ; dashed-curves, ECS results [9, 50]; dash-dotted curves, CCC 
results [26, 48]. The experimental results are the absolute measured values of Roder et 
al [47,48]. Present results have been multiplied by a factor 0.5 (see text). 
Figure 1(c). TDCS in coplanar equal-energy-sharing with two electrons emerging on 
opposite sides of the direction of the incident electron with equal angle 9a and energy 
Ei = 15.6 eV . Continuous curves, present results ; dashed-curves, ECS results [9, 50]; 
dash-dotted curves, CCC results [26,48]. The experimental results are the absolute mea- 
sured values of Roder et al [47, 48]. Present results have been multiplied by a factor 0.5 
(see text). 

Figure 2(a). Same as in figure 1(a) but for 17.6 eV incident electron energy. The exper- 
imental results are the recent absolute measured values of Roder et a/ [51] and the CCC 
results are as in [29]. Here the present results arc free from any multiplicative factor. 
Figure 2(b). Same as in figure 1(b) but for 17.6 eV incident electron energy. The 
experimental results are the recent absolute measured values of Roder et al[51] and the 
CCC results are as in [29] . Here the present results are free from any multiplicative 
factor. 

Figure 2(c). Same as in figure 1(c) but for 17.6 eV incident electron energy. The 
experimental results are the recent absolute measured values of R5der et a/ [51] and the 
CCC results are as in [29]. Here the present results are free from any multiplicative factor. 

Figure 3. Same as in figure 2(a) but for 19.6 eV incident electron energy. The relative 
measured results of Roder et al [47-49] are normalized as in [9] . 

Figure 4. Same as in figure 3 but for 25 eV incident electron energy. The CCC results 
are from [28, 49] 

Figure 5. Same as in figure 3 but for 30 eV incident electron energy. 
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